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ABSTRACT 



Aims. We investigate the possibility of using the X-ray telescope (XRT) on board the Swift satellite to improve the current accuracy 
of the intra-cluster medium (ICM) temperature measurements in the region close to the virial radius of nearby clusters. 
Methods. We present the spectral analysis of the Swift XRT observations of 6 galaxy clusters and their temperature profiles in 
the regions within 0.2-0.6 r 2 oo- Four of them are nearby famous and very well studied objects (Coma, Abell 1795, Abell 2029 and 
PKS0745-19). The remaining two, SWJ1557+35 and SWJ0847+ 13, at redshift z=0.16 and z=0.36, were serendipitously observed by 
Swift-XRT. We accurately quantify the temperature uncertainties, with particular focus on the impact of the background scatter (both 
instrumental and cosmic). We extrapolate these results and simulate a deep observation of the external region of Abell 1795 which is 
assumed here as a case study. In particular we calculate the expected uncertainties in the temperature measurement as far as r 2 oo- 
Results. We find that, with a fairly deep observation (300 ks), the Swift XRT would be able to measure the ICM temperature profiles 
in the external regions as far as the virial radius, significantly improving the best accuracy among the previous measurements. This 
can be achieved thanks to the unprecedented combination of good PSF over the full field of view and very accurate control of the 
instrumental background. 

Conclusions. Somehow unexpectedly we conclude that, among currently operating telescope, the Swift-XRT is the only potentially 
able to improve the current accuracy in plasma temperature measurement at the edges of the cluster potential. This will be true until 
a new generation of low-background and large field of view telescopes, aimed to the study of galaxy clusters, will operate. These 
observations would be of great importance in developing the observing strategy for such missions. 



1. Introduction 

Galaxy clusters form by the hierarchical accretion of cos- 
mic matter. They reach the virial equilibrium over a volume 
that defines the regions where the pristine gas accretes on the 
dark matter (DM) halo through gravitational collapse and is 
heated up to millions degrees through adiabatic compression 
and shocks. The end products of this accretion process ex- 
hibit i n the X-ray band sim i lar radial profiles of surface bright- 
ness dVikhlinin et alJ Il999t INeumannl l2005b lEttori & Balestral 
2009), plasma temperature (Allen et alJ I200U IVikhlinin et alJ 
I2Q05I: ILeccardi & Molendil I2008 
distribution ( Pointecouteau et ail 



and gravitational mass 
2005, e.g.). The measurement 



of the properties of the ICM have been enormously im- 
proved thanks to the arcsec resolution and large collecting 
area of Chandra and XMM-Newton, but still remain possible 
only where the X-ray emission can be well resolved against 
the background (both instrumental and cosmic). While the 
X-ray surface brightness an d gas density can be estimated 
in few cases above 0.7r2oo (IVikhlinin et all 1 1999t INeumannl 
I2005t lEttori & Balestral l2009l e.g.), the ICM temperature, re- 
quiring more than an order of magnitude in net counts than 
the surface brightness to be firmly measured, can be reason- 
ably well con strained up to a fraction (~ 0.5 - 0.6) of the 
virial radius (IVikhlinin et alJl2006t ILeccardi & Molendil l2008t 
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lEttori & Molendi 201(1 lArnaud et al.l [20 1 OL e.g.). The regions 
not-yet observed are expected to retain most of the infor- 
mation on the processes that characterize the accretion and 
evolution within the cluster of the main b aryonic component 
dRoncarelli et al.ll2006t Ikasheed et al. 2010). It is therefore cru- 
cial to obtain direct measurements of the cluster properties at 
these large radii where very important processes for the evolu- 
tion of the clusters take place. Very recently, Suzaku, thanks to 
its low background and high sensitivity, has been able to map 
roughly (i.e. with a spatial resolution limited to >4') the regions 
close to the virial radius, providing the firs t estimate of the gas 
temperature in 6 objects (Fuiita et al. 2008; George et al. 2009; 
Reip rich et al1 | 2009[ iBautz et alJl2009t iKawaharada et alJboiot 
iHoshino et al.ll2010» . T he aim of this paper is to show that the 
X-ray tel escope (XRT) dBurrows et al.ll2005h on board the Swift 
satellite dGehrels et alJ [2004) can improve the accuracy of these 
measurements with fairly deep observations. To this end, we 
present the first spectral analysis of the archived Swift observa- 
tions of nearby clusters (Section|2]i. Since this is a non-standard 
analysis for the Swift - XRT data and is used and presented here 
for the first time, we discuss in detail the techniques adopted. 
In particular in Section [3] we describe the procedure we devel- 
oped to estimate the background and its systematic uncertainty. 
To calculate how much this uncertainty affects the temperature 
measurement we performed a series of Monte Carlo simulations 
of thermal spectrum at different level of surface brightness (de- 
scribed in Section|4]). Finally in Section|5]we used these results 
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Fig. 1. Images for the six objects in our sample, with detected sources excised. Data are smoothed with a 30"Gaussian filter. 
Superimposed are the circular extraction regions used for the spectral analysis. Circular regions around the detected sources are 
excised. 

Table 1. Cluster sample. Mean spectroscopic temperatures and redshift of the first four clusters are from the references reported in 
the fifth column, through a query to BAX archive (http : //bax . ast . obs-mip . fr /); the remaining 2 are calculated from our data 
(see text). V200 are estimated by adopting the scaling relations dArnaud et alj|2005h and here are expressed in units of arcmin and 
Mpc. r max is the center of the most external annulus for which we measured the temperature, here expressed in units of 1-200 ■ 



Name 


RA, Dec 


z 


kT 


references 


r 200 


"max 


Exp. 


N H 




[deg] 




[keV] 




[', Mpc] 


[ r 20o] 


[ks] 


10 2( W 2 


Coma 


194.9392 +27.9429 


0.023 


8.25+0.01 


Arnaudetal. (2001) 


78.2 [2.18] 


0.17 


43.7 


0.80 


Abell 1795 


207.2183 +26.5903 


0.062 


6.12+0.05 


Vikhlinin et al. (2006) 


25.7 [1.84] 


0.41 


20.3 


1.32 


Abell 2029 


227.7336 +5.74440 


0.077 


8.47+0.09 


Vikhlinin et al. (2006) 


26.7 [2.16] 


0.43 


42.3 


3.25 


PKS0745-19 


116.8799 -19.2948 


0.102 


7.97+0.28 


Arnaud et al. (2005) 


18.1 [2.06] 


0.58 


63.1 


41.8 


SWJ1557+3530 


239.4287 +35.5073 


0.153+0.006 


6.79+0.25 


this work 


11.5 [1.89] 


0.62 


180.1 


0.20 


SWJ0847+1331 


131.9550 +13.5278 


0.358+0.005 


6.02+0.34 


this work 


5.3 [1.56] 


0.50 


203.1 


3.23 



to properly and robustly evaluate the expected uncertainties on 
the gas temperature measurements at r2oo El on a simulated 300 
ks observation of the well studied cluster Abell 1795. 

Throughout this paper we assume Ho=70 km s _1 and Q,t=0.73 
and O m =0.27, which are the default values in the XSPEC(vl2.5) 
software. All errors are quoted at 68% confidence level for one 
parameter of interest, unless otherwise specified. 



1 The radius that defines the sphere enclosing a mean cluster density 
that is 200 times the critical value at the cluster's redshift. Here we use 
r 2 oo and virial radius indifferently and we calculate r 2 oo using the scaling 

relations given by I Arnaud et al] J2005h r 20 o = 1714^2— [kpc] where 

hz= V((l +z) 3 Sl,„ + n A ) 



2. Data reduction and analysis procedures 

2.1. The sample 

The X-ray telescop e (XRT) on board the Swift satellite 
dGehrels et al.ll2004h . uses a Wol ter I mirror set, orig inally de- 
signed for the JET-X telescope (ICitterio et al.l fl994). to focus 
X-rays ( 2- 10 keV) onto a X MM-Newton/EPIC MOS CCD 
detector (Burrows et al. 2005). The effective area of the tele- 
scope (~ 150 cm 2 at 1.5 keV) is ~ 3. smaller than one sin- 
gle XMM-Newton MOS. The PSF, similar to XMM-Newton, 
is ch aracterized by an ha lf energy width (HEW) of ~ 18"at 1.5 
keV dMoretti et al.l2 005). As part of the Swift scientific payload, 
since the start of the mission (November 2004), Swift-XRT has 
been mostly used to observe GRB afterglows and other variable 
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Fig. 2. Upper panels: The projected temperature profiles of the six clusters o f our sample. For the first four objects, Swift XRT 
measurements (black points) are com pared with XMM-Newton (blue points) ISnowden et alj d2008l) . In the Coma plot we also 
report the BeppoSA X numbers from l De Grandi & Molendil(l2002l) . In the Abell 2029 and Abe ll 1795 plots we rep ort the Chandra 
measurements from IVikhlinin et a l. (2006). In the PKS0745-19 we report the Suzaku points dGeorge et al.ll2009l) . In the case of 
SWJ1557+35 and SWJ0847+13 we compare measurements obtained with locally evaluated and modeled background. M iddle 
panels: Fluxes in the 0.3-10. keV band in cgs deg~ 2 units compared with XMM-Newton values from lSnowden et al] (120081) . in the 
first four cases, and with the value after local background subtraction in the remaining two. Lower panels: For each extraction 
region we plot the relative value of the 0.5-2.0 keV band flux for the cluster (black), XRB (green) and NXB (red). 
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sources. Galaxy clusters are observed mostly for calibration pur- 
poses or serendipitously. 

We cross-correlated the XRT archive updated on April 2010 with 
the BAX database (http://bax.ast.obs-mip.fr/) search- 
ing for cluster observed with more than 20,000 events registered 
in the 0.7-7.0 keV energy band. The sample consists of four fa- 
mous and very well studied objects (FigQ]Tab.[T]). For these clus- 
ters the XRT archive data allowed us to measure the projected 
temperature profile as far as the edge of the telescope field of 
view corresponding to 0.2-0.6 r2oo- 

We added 2 more clusters to this sample; these are 
SWJ1557+3530 and SWJ0847+1331 at redshift z=0.15 and 
z=0.36 respectively (according to X-ray spectral analy- 
sis, see below ), which were serendipitously observed in 
the field of the follow-up observations of GRB090409 
and GRB051016B and are t he two highest si gnal ob- 
jects of the SXCS catalog (IMoretti et all 120071 Moretti 
et al. in preparation). They are already catalo ged as 



MaxBCGJ239.42665+35.50827 dKoester et all l2Q07h (pho- 
toz=0.155) and as WHL J084749 3+133140 dWen et"al1 l2009) 
(z=0.36) respectively. They are X-ray bright (Lx ~ 3 and 7 10 44 
erg s~l in the 0.5-2.0 keV rest frame energy band) with av- 
erage temperatures of 6.0 and 6.8 keV respectively. They can 
be considered rich also form the point of view of the optical 
classification having 35 and 42 galaxies associated according 
the MaxBCG and the WHL catalogs respecitvely. The cataloged 
redshift measurements are in very good agreement with ours (see 
Tab. [!)■ Differently from the first four, these two objects do not 
entirely fill the field of view allowing us to perform a consistency 
check of our background model through a comparison with the 
local background (see below). 



2.2. Data Reduction 

Data reduction was performed using the standard software 
(HEADAS software, v6.8, CALDB version 20091130, Nov 
2009) and following the procedures described in the instrument 
user guide Q. At variance with these, we excluded the external 
(Detx > 90 and Detx <510 ) CCD columns which are affected 
by the pres ence of out-of-time -events from corner calibration 
sources (see iMoretti et al. 2009, for a detailed map of the CCD 
and a discussion on the XRT background). This left us with a 
nominal field of view of 16.5' x 18.9' (0.087 deg 2 ). Different ob- 
servations of the same objects and relative exposure maps were 
merged by means of the the extractor and f arith tasks of the 
HEADAS software respectively. 

Before performing the spectral analysis, we ran the CIA0 
wavedetect and eliminated all the events within the circles cen- 
tered on the detected source positions with radius such that the 
PSF surface brightness equals the background (typically < 10 
pixels < 23"). Vignetted and non-vignetted exposure maps were 
built consequently accounting for the excised regions. 



2.3. Spectral Analysis 

We measured the temperature profiles, binning the data in pro- 
jected annuli of amplitudel.0'-2.5' moving from central to the 
external part of the image (Fig [TJ. The statistics of the extracted 
spectra strongly vary going from the center, where the extraction 
region size is determined only by the PSF toward the outskirts of 
the sources, where the size of the regions is determined by ensur- 
ing a minimum of 200 source counts in the 0.7-2.0 keV energy 




http://heasarc.nasa.gov/docs/swift/analysis/documentation 



2 4 6 8 10 12 

XMM temperature [keV] 

Fig. 3. Swift-XRT and XMM-Newton temperature measurement 
comparison. The continuous line is the linear fit to the data, 
while the dotted line is the expected distribution. In the inset the 
scatter of the measurement is shown with the thick black line. 
The best fit with its uncertainties is shown with thin black lines 
and grey area, while the expected scatter is drawn with red dotted 
line. 



band. Because the XRT standard software does not include any 
specific task for the spectral analysis of the extended sources, 
we developed a procedure made by combination of ad-hoc IDL 
routines and Ftools scripts to calculate the vignetting and the 
exposure corrections and to estimate the background. 
For the vignetting calculation, we, first, created a vignetting map 
at 1.5 keV, dividing the exposure map by the un-vignetted expo- 
sure. Then, we created the maps at differe nt energies, u s ing th e 
vignetting analytical description given by Moretti et al. (2004). 
Finally we used the exposure and vignetting maps weighted by 
the cluster surface brightness profile to calculate the effective 
time and the vignetting correction, for each annulus. As the vi- 
gnetting depends on energy, in order to take into account for 
these corrections in the spectral analysis we opportunely modi- 
fied the nominal ARF file. 

The evaluation of the background is surely the most important 
and delicate step in the procedure. This is primarily because the 
virial regions (R200 > 15') of nearby massive clusters subtend 
the entire field of view. Therefore, the background cannot be 
estimated locally. For our purposes we consider the total XRT 
background as the sum of an instrument (NXB) and a cosmic 
(XRB) component. The latter one is, in turn, the sum of a galac- 
tic (GXRB) and extragalactic (CXRB) component, with very dif- 
ferent spectral and spatial characteristics. 

We evaluated the instrumental background (NXB) using the data 
collected, during the observations, in the regions of the detec- 
tor which are not exposed to the sky (NESR, Not-Exposed-Sky- 
Regions). These are four different small regions (2507 pixels 
each) close to the CCD bou ndary and delimited by the field of 
view and the corner sources (IMoretti et al. 2009). 
To estimate the XRB we used a statistical approach: we made 
use of a large number (135) of deep GRB follow-up observa- 
tions to study the XRB characteristics. Because GRB are uni- 
formly distributed, these are deep exposures on random positions 
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of the sky, totally uncorrelated with already known bright X-ray 
sources, and provide us with a data set which is very well suited 
to statistically characterize the XRB. In the following we will 
refer to this dat a set as the blank fi elds (BFs). A similar data-set 
was used by iMoretti et alj J2009) to perform an absolute mea- 
surement of the CXRB in the 1.5-7 keV range. Here, we found 
that all the BF spectra, in the 0.7-7.0 keV energy band can be 
well described by the sum of a bremsstrahlung and an absorbed 
power law, representing the GXRB and CXRB respectively, 
We defer the detailed description of the instrumental and cos- 
mic background analysis together with the full evaluation of the 
systematic errors in temperature measurements to the next two 
sections (Sect. [3] |U, while here we focus on the remaining de- 
tails of the analysis and the results. 

As high energy limit we use 7.0 keV, because, between 7 and 
10 keV, the instrumental background is significantly higher due 
to the presence of Nickel and Gold lines. We set the low en- 
ergy limit at 0.7 keV both because around 0.5 keV response 
matrix calibrations are more problematic due t o the presence 
of uncorrected charge traps dGodet et~a"fl l2008h . and because 
we found that, below 0.7 keV, the XRB is less reproducible 
due to the presence of s ome extra local components (see also 
iKuntz & Snowder] (I2000I) '). In the case of Coma, both brightness 
and temperature are so high that we easily extended our analy- 
sis up to 10.0 keV. In the case of PKS0745-19, which is on the 
Galactic plane, we used the 1.0-7.0 keV band to reduce the un- 
certainty on the background estimate which is higher at lower 
Galactic latitudes. 

We used XSPEC(vl2.5) Cash statistics to fit the clus- 
ter+XRB spectra from the annular regions with a 
bremss+wabsx(pow+apec) model, where the absorbed 
APEC model represents the cluster emission, using the NESR 
spectra as background. We used literature red-shift values 
(http://bax.ast.obs-mip.fr/) for all the clusters but 
for SWJ1557+3530 and SWJ0847+1331 which are not yet 
cataloged. In these cases we measured the redshift using the 
Fe-K line in the X-ray spectrum from the central 1' radius. 
We left the metallicity parameter free to vary only when the 
estimated counts of the sources exceed 5,000, freezing it to the 
0.3 Z in the other cases. Concerning the Galactic absorption, 
for each cluster, we took advantage of the good statistics in the 
three central annuli, and we let the N# absorbing column vary 
in the fit. As we found, in all the cases, a <2<x consistency of 
the b est fit value with th e value derived from the HI Galaxy 
map ( Kal berla et alJl2005l) . we froze the Galactic N# absorbing 
columns to these values. In three cases these estimate coincide 
with the XMM bes t fit values (Com a, Abell 1795, PKS0745-19 
differences <3%) dSnowden et al.l l2008h . while, in the case of 
Abell 2029, the difference is -20%. 

We tried to fit the cluster data together with XRB, but the statis- 
tics of the data did not allow us to constrain the XRB parameters 
together with the cluster ones, even when the XRB is a signifi- 
cant fraction of the total signal. Instead, for all our spectra, we 
froze the four XRB parameters to the median values of the BF 
distributions, just re-normalizing them by the area of the extrac- 
tion region. We calculated the systematic uncertainties of this 
approach reproducing the same procedure on a large sample of 
simulated spectra (Sect. [4]). Finally we calculated the total error 
as the quadratic sum of the statistical error (the one provided 
by XSPEC) and the systematic one. The results are shown in the 
Fig. [2] and reported in details in Tab. [2] 

We note that in the case of the PKS0745-19 the systematic er- 
rors are probably underestimated, due to the fact that this cluster 
lies on the Galactic plane, whereas our XRB statistical analy- 
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Fig. 4. The Aitoff projection of the Galactic coordinates of the 
non resolved emission in the blank fields (BFs). Upper Panel: 
the CXRB flux in the 2.0-10.0 keV band. Lower Panel: the 
GXRB flux in the 0.5-2.0 keV band. 



sis is suitable for the extragalactic sky (Galactic latitude > 20°). 
We mitigated this problem limiting the spectral analysis for this 
particular cluster to energies higher than 1.0 keV. 
In Fig. [3] we plot the XRT temperature measurements of the 
first four clusters, c ompared with XMM-Newton results from 
ISnowden et al.l(l2008l) : we generally found a good agreement be- 
tween the two instruments without any appreciable bias or sys- 
tematic trend as shown by the linear fit. On the other hand, calcu- 
lating the scatter as the variance of the (XRT)-(XMM-Newton) 
difference distribution, divided by the errors, we measured a 
scatter which is slightly larger than the expected: 1.32+0.17 in- 
stead of 1.0. Most of this scatter can be ascribed to the outer 
regions of the cluster PKS0745-19, for which, as said, the errors 
are underestimated. Indeed excluding them lowers the scatter to 
1.20+0.18, consistent with the expectation at 1.1 <x level. We 
obtained the same result (1.22+0.20) artificially increasing the 
PKS0745-19 errors by a factor 1.4. In this comparison we per- 
form ed the spectral anal ysis exactly in the same regions used 
by (ISnowden et al.l l2008). while data plotted in Fig. [2] and listed 
in Tab. [2] are in (slightly) different regions, modified in order to 
optimize XRT data statistics. 

Due to small apparent size and the good statistics 
SWJ1557+3530 and SWJ0847+1331 are the best cases, in 
the XRT archive, to perform a consistency check between 
temperature measurement performed with a locally evaluated 
and a background model: results are in good agreement, as 
shown in the last two panels of Fig. [2] 
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Table 2. Extraction regions and spectral analysis results. For each annulus we report the radii, the area, the effective exposure time 
accounting for exposure maps and vignetting, the mean cluster surface brightness, the total events, the source/total ratio and the best 
fits values. The statistical error contribution to the total temperature errors are reported in parenthesis. 
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3. Background Analysis 

To study both the NXB and XRB we made use of a large sample 
of blank fields (BFs). We selected the 135 (31 at low Galactic lat- 
itude) GRB follow up observations from January 2006 to April 



2009 with a nominal exposure time longer than 10 ks and shorter 
than 300 ks (Fig. [4}. We reduced these data and eliminated de- 
tected sources, following the same procedure we used for cluster 
data and described in Sect. 12.31 
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Fig. 5. Upper panel: The raw spectrum (P.I. channel distribu- 
tion; 1 PI~ 10 eV ) of a large number of blank fields (BF, in 
grey) compared with the spectrum registered in the NESR re- 
gions in black. Middle panel: The ratio NESR/BF gives the rel- 
ative importance of the NXB as function of energy in absence 
of bright sources. Bottom panel: The 4 NESR spectra show not 
negligible differences. 



3. 1 . Instrumental background 

In the Swift-XRT CCD, we have the opportunity of estimating 
the NXB directly from the signal registered in the NESR during 
the observations. This method allows to optimally map the NXB 
time variations which occur due to changes in the satellite en- 
vironment at any time scales from minutes to years and that are 
very difficult to be accounted for otherwise. Indeed, as shown 
in the Fig. [5] the raw spectrum registered in the NESR regions, 
once re-normalized for the area, accurately reproduces the sig- 
nal registered in whole the field of view beyond 7 keV, where the 
instrument background completely dominates the BF signal. On 
the contrary NESR are in the CCD corner s and do not map th e 
NXB spatial variations within the detector. iMoretti et all d2009). 
using the data collected during the unique (and unrepeatable) ob- 
servation performed with the camera shutter closed (September 
2007), found a linear gradient in the NXB signal along the ver- 
tical direction of the CCD, with the bottom regions being 20% 
fainter than the top regions. This gradient was found achromatic 
within the statistical errors (Icr ~5%). Here, we used all the BFs 
stacked data (6.3 Ms in total) to study the NXB spatial pattern. 
First, comparing the 4 NESR spectra we found that while the 
upper regions (NESR 1,2) and lower (NESR 3,4) are consistent 
between each other both in intensity and spectral shape, the up- 
per ones are significantly different from the lower ones (Fig. [5]). 
Then, to study the NXB spatial pattern we used the unresolved 
signal from all the BFs stacked data in the 7.0-10. keV energy 
band. As said, in absence of bright sources, in this range, the sig- 
nal is almost purely instrumental, dominated by the Nickel and 
Gold fluorescence lines (Fig. [5]). From these data we confirmed 
that the NXB has a linear vertical gradient with slope very well 
consistent with the one found by More tti et alj (|2009). The scat- 
ter in the stacked BF 7-10. keV image measured in 50x50 pixel 
cells is ~ 8%, (1.5% statistical); correcting the image by the lin- 
ear gradient reduces this scatter to ~ 2% (Fig. [6]). If we assume 
that the NXB linear gradient is achromatic this means that, start- 
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Fig. 6. Upper-left: The CCD map in the 7.0-10.0 keV energy 
band from 6.3 Ms of BFs, normalized to the central raw to show 
the NXB vertical gradient; Lower-left: the residual distribution 
calculated in 50x50 pixel cells. Upper-right: The CCD map in 
the 7.0- 10.0 keV energy band from 6.3 Ms of BFs, normalized to 
the central raw to show the NXB vertical gradient, corrected by 
the linear approximation. Lower-left: the residual distribution 
calculated in 50x50 pixel cells with the linear gradient correc- 
tion applied. 



ing from the NESR signal, we can recover the NXB with ~ 2.3% 
accuracy (1 cr, plus the statistical uncertainty) in any position of 
the detector. The assumption of achromaticity is justified by the 
analysis results of the shutter-closed observation. 

3.2. Cosmic background 

In order to study the statistical properties of the cosmic com- 
ponents of the background we performed the spectral analysis 
of the 135 BFs unresolved signal, using the NESR spectrum 
as instrumental background. Freezing the value of the Galactic 
absorbing Hydrogen column, for each field, to iKalberla et al.l 
(2005) values, we found that the unresolved emission in the 0.7- 
7.0 keV can always be well modeled by a thermal component 
(bremsstrahlung) plus an absorbed power law, accounting for 

GXR B and CXRB re s pectiv ely. 

From ISnowden etaD (1 19981) and dKuntz & Snowdenll2000l) we 
know that a physically motivated emission model of the XRB be- 
low 1 keV should consist of, at least, two components, the first 
of Galactic origin, the second from local hot bubble. To avoid 
this complexity in the modeling of the data, we limited our anal- 
ysis to energies higher than 0.7 keV. Although not completely 
physically motivated, using a simple Bremsstrahlung provided 
us an accurate phenomenological description of the data in this 
energy range with a limited number of parameters. This is suit- 
able to our goal, whereas a physical description of the galactic 
thermal emission is beyond the scope of the present work (and 
the quality of our data at low energies). 
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Fig. 7. Upper panel: The reduced^ 2 distribution of the 135 fits, 
here split in high (> 20 °) and low Galactic latitude fields. Lower 
panel: The 104 XRB models (grey lines) with the median spec- 
trum in evidence (thick black line). For this we over plot the two 
components: the Galactic bremsstrahlung (red dashed ) and the 
extragalactic power-law (blue dashed) 



Table 3. The median of the BF spectrum parameter distribution 
together with the 16th and 84th percentile values. 





median 16th % 84th % 


B. norm [phot s 'cm 2 deg 2 kev '] 
B. kT [keV] 

P.L. norm [phot s _1 cm _2 deg _2 kev _1 ] 
P.L. slope [] 


0.023 0.0075 0.069 
0.27 0.21 0.42 

0.0021 0.0018 0.0024 
1.20 1.08 1.34 
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Indeed we found good x 2 values for all the 135 fields (upper 
panel of Fig.[7j. Splitting the sample in low (< 20) and high lati- 
tude Galactic fields we found that fitting the latter always yielded 
X 2 <1.3, while the formers present slightly larger scatter. The 
soft thermal component typically contributes to ~50% of the to- 
tal emission <1 keV, while it is negligible beyond 1.5 keV. On 
the other hand, as we will see in the next Section, neglecting 
it, could significantly affect the temperature measurement in a 
regime where the cluster emission is comparable. As expected, 
we found that high and low latitude fields have different statis- 
tical properties, the latter sample presenting a larger scatter in 
the parameter distributions. As the cluster in the present sample 
are all (but PKS0745-19) observed at high Galactic latitude we 
restricted our analysis to the extra Galactic fields. The parameter 
space of the best fit values of the 104 high Galactic latitude BFs 
is shown in Fig. [8] We used these data in two ways. First, from 
these distributions we derived the median spectrum (lower panel 
°f Fig-El Fig.[8]and Tab. [3]) used as XRB background component 
in the cluster spectral analysis (see Section 12.3b . Then, we used 
these models as input templates for the simulations employed to 
quantify the systematic temperature uncertainties (Section|4]i. 
Approximately ~20% of the unresolved emission, in the XRT 
images, is due to stray li ght contamination f rom XRB sources 
outside the field of view (Morettietal3[2009). This contamina- 
tion does not significantly affect our procedures, but it should be 



Fig. 8. Upper panel: The bremsstrahlung parameter space of 
the 104 high Galactic latitude fields. The red dotted lines repre- 
sent the 16th and 84th percentile of the distributions and the red 
triangle the median. The red triangle represent the median values 
adopted as XRB model (Tab. Lower panel: The power-law 
parameter space, with same notations. 



taken into account to compare the normalization of our BF spec- 
tra with the absolute measurements of the CXRB and GXRB. 



3.3. Time variability 

Time variations of both cosmic and instrument background level 
are observed in the X-ray telescopes. In the Swift XRT par- 
ticle flux variations in the satellite environment affecting the 
NXB, are monitored and accounted for by the NESR data (see 
Sect l3.1b . More insidious are the variations due to the solar wind 
exchange (SWCX) producing a diffuse photon emission at C, O, 
Mg, Ne energies with different timescales from seconds to days. 
We quantified the impact of SWCX on Swift XRT observations 
studying the light curves in different energy bands both for the 
cluster data and for the BFs. 
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Fig. 9. The light (upper-left panel) and color (lower-left) curve 
of the unresolved signal in the SWJ1557+35 observations. In the 
right panels we plot the corresponding distribution together with 
the gaussian fit. The observed scatter can be completely ascribed 
to the statistical error of the measurements. 

We studied the light curve in the soft band (0.7-2.0 keV) and 
the color curve (0.7-2.0/0.7-7.0 keV) of the source free regions 
of the 2 serendipitous cluster observations (SWJ1557+35 and 
SWJ0847+13). We excluded the sources as we did for BFs and 
binned the data in ~ 1,500 s time intervals 0. These choices en- 
sured the necessary statistics to study the flux and spectral vari- 
ations, as the signal registered in a BF is ~ 0.1 count per sec- 
ond in the 0.7-7.0 keV energy band over the whole detector. 
Results relative to SWJ1557+35 are shown in Fig. [9] The to- 
tal exposure time of the observation is 200 ks which is split in 
149 segments. We considered only the 1 15 intervals with a dura- 
tion longer than 1,000 s covering more than the 90% of the total 
duration (185 ks). We found that both the flux and the color scat- 
ter are well consistent with the statistical one. In other words, in 
this data set, we did not find any detectable signature of SWCX. 
For SWJ0847+13 we found similar results, while for the four 
remaining objects the results of this analysis are not conclusive, 
due to the very high level of the cluster signal over the entire 
field of view. 

We extended this kind of analysis to all the BFs. To get rid 
of the cosmic variance, we normalized each observation to its 
median value and compared the residuals. We found that only ~ 
1.5% of all the time intervals considered can be identified out- 
liers Fig. [10] Given these findings, we decided to neglect the 
effect of the SWCX in both our spectral analysis and in our sim- 
ulations. 

4. Systematic uncertainties in temperature 
measurement 

In order to evaluate the systematics in the temperature measure- 
ment we used XSPEC(vl2.5) to simulate a set of cluster+XRB 
bremss+wabsx(pow+apec) spectra, as observed by the Swift 
XRT 

To simulate the XRB we used BFs, but we had to consider two 
further complications. First, CXRB normalization and spectral 

3 Due to the low orbit of the satellite, the Swift XRT observations are 
composed by segments of ~ 1,500 seconds for each orbit (5700 s.). 
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Fig. 10. The residual distribution of the flux and color XRB 
curve of a collection of BFs. We put in evidence the expected 
Gaussian distribution from statistical uncertainties and the de- 
partures from this. We found that only ~ 1.5% of all the time 
intervals considered can be considered outliers. 



slope are expected to vary according the image flux limit: the 
spectrum of the unresolved CXRB emission in deeper images is 
expected to be fainter and harder (Mo retti et al.ll2003h . Second, 
the CXR B variance depends on the size of the extraction re- 
gions dKushino et alJl2002l:lRevnivtsev et~aT1l2008l:lMoretti et al.1 
2009). To account for these effects we split the BF sample in 
five exposure time bins (10-30, 30-50, 50-75, 75-125, 125-200 
ks) and for each BF we extracted the XRB spectrum from re- 
gions of different sizes (1.0, 0.5, 0.25 and 0.1 times the entire 
field of view). 

To simulate the cluster emission we used 37 different cluster 
temperatures (in the range 1-10 keV, step 0.25 keV) with the 
same flux (8xl0~ 13 erg s _1 crrT 2 ) and exposure time (100 ks) in 
order to collect a number of source photons to make the statis- 
tical errors negligible (> 10000 source counts ). Then, for each 
spectrum we modified the source/background ratio, varying the 
BACKSCAL keyword, in order to simulate different values of 
surface brightnesses: 10 steps between 2xl0~ 12 and 8xl0" n 
erg s crrT 2 deg~ 2 ). For each step of this 37x10 grid we sim- 
ulated 200 different realizations letting the metallicity and Nh 
normally and randomly varying around the mean values of 0.3 
Z and 3xl0 20 cirT 2 respectively with a scatter of 0.06 (20%) 
and 6xl0 19 (20%). Thus, for each simulated cluster we summed 
an XRB model randomly choosing one of the BF model cycling 
over the exposures x region size grid: this resulted in a total of 
200 spectra x (37x10) clusters (x 5x4) XRB grid. Finally, to 
each simulation, we added the instrumental background stack- 
ing a collection of the NESR data, randomly re-normalizing it 
by a gaussian deviation of 2.3% to account for the uncertainties 
in the reproduction of the spatial pattern (Section [3~Tl i. We used 
all the BF observed in 2009, for a total of ~1.5 Msec, in order 
to maximize the statistics and, at the same time, to avoid parti- 
cle background vari ation on ~ 1 year time scale , as observed in 
iMoretti et alj (|2009), which are probably due to variation in the 
Solar activity. 

We fit the simulated data using the same procedure we 
used for our real data as described in Sect 12.31 We used a 
bremss+wabsx(pow+apec) model, freezing the XRB param- 
eters to the BF XRB median spectrum (Sect. [3~2l ). just normal- 
ized for the input area. Moreover, we froze the Nh and metal- 
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licity to the mean values (0.3 Z and 3xl0 2() cirT 2 respectively 
) and we used the GRB 090618 follow-up observation NXB as 
background. 

As an example of the simulation outputs, in Fig. Q~T] we plot 
the distribution of the results of the fit on the 200 realizations 
of APEC spectra with different input temperatures T,, for two 
different values of surface brightness (high surface brightness in 
the higher panel, low surface brightness in the middle panel). In 
this case the background was extracted from BF with exposure 
time within 75-125 ks and size ~ 0.5 the field of view. It is im- 
mediately evident that at lower values of SB the scatter of the 
simulation outputs is significantly higher. We used these outputs 
to calculate the systematic uncertainties in our cluster tempera- 
ture profiles in the following way. We did not calculate the error 
of a measurement T,„ simply as the scatter of the results of the fit 
of the 200 realizations at temperature T m . Instead we calculated 
it from the distribution of the "true" temperatures T, that have 
a non null probability of being measured as T m . When the SB 
is high (> le-11 erg s _1 cirT 2 deg~ 2 ) the systematic scatter of 
the measurement T, is very small and the probability that a true 
temperature T, » T„, would yield a T m measurement is very 
low. For example in the case shown in the upper panel of Fig.fTTI 
with a SB ~ 5 10~ n erg s _1 crrT 2 deg~ 2 the probability that a 
cluster region with a true temperature T, >4 would be measured 
3 keV is almost null. On the contrary, when the SB is low and 
the background is comparable or higher than the cluster signal, 
the scatter of the measurement can be very large, especially at 
high temperatures. In the middle panel of the same figure, we 
show the same example for a SB which is a factor 10 lower (5 
10~ n erg s _I cirT 2 deg~ 2 ): in this case the probability that T; >4 
would be measured 3 keV is not negligible. Thus, to calculate the 
T„, error, for each T, of the simulation grid, we calculated the T, 
output distribution value in T,„; then we calculated the errors as 
the 16th and 84th percentile of this distribution (lower panel of 
FigHB. 

The systematic errors reported in Tab. [2] and plotted in Fig. [2] 
are calculated in this way. In most of the spectra of the present 
sample, the systematic errors do not strongly affect the total er- 
ror budget of our measurements as the surface brightness of the 
observed regions is higher than 10~ n erg s _1 cirT 2 deg~ 2 . The 
external regions of the two more distant clusters are the only 
cases where SB is lower; but here the statistical errors are very 
large too. However the analysis of the temperature measurement 
systematic errors is mandatory to realistically calculate the ex- 
pected errors in the outer regions assuming to have a good statis- 
tics data-set as in our simula tion of Abelll795. 
ILeccardi & Molendl d2007l) shown that, using the most com- 
mon likelihood estimators, the measurements of the tempera- 
tures from outer regions of galaxy clusters are strongly biased. 
We note that our way of calculating the temperature measure- 
ment uncertainties represents a refinement of the usual procedure 
used in the literature. Indeed, the current approach is a particular 
case of our procedure, given by the simplifying assumption that 
the measured value T,„ coincide with the true value T,. In the re- 
maining we will refer to this way of calculating the error as HI, 
while we will call our refined procedure H2. 



5. A1795: a case study 

Chandra and XMM-Newton deep observations showed that 
Abell 1795 is a dy namically relaxed system within rsoo 
(IVikhlinin et al.ll2006h . altho ugh a 50 kpc (40 ") long X-ray fil- 
ament in the core is present (Fabi an et alj|2001[) . Thanks to the 
lower background, Suzaku, for the first time, succeeded in ex- 
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Fig. 11. Upper panel: The simulation outputs relative to an 
APEC spectrum with kT=l-7 keV 100 ks observation extracted 
from a region size of 0.5 the field of view in a high signal to back- 
ground ratio regime (SB = 5 10~ n erg s _1 cirT 2 deg~ 2 ). Different 
colors indicate the distribution of the results of the fits of differ- 
ent input temperatures (200 realizations each). Middle panel: 
The same of the upper panel in the case of low signal to back- 
ground ratio (SB = 5 10~ 12 erg s" 1 cm -2 deg~ 2 ). Lower panel: 
The distribution of the T, scatter values measured at T m =3 keV 
relative to the two different values of surface brightness; the in- 
tervals between 16th and 84th percentile of both distributions are 
in evidence. 



tending the observations in the radial range 1-500 and r2oo finding a 
significant asymmetry between the northern and southern part of 
the cluster. Departure form the hydrostatic equilibrium, with the 
plas ma in-falling from the northern skirt is a possible explana- 
tion dBautz et al.l2009l) . Suzaku observation sensitivity is limited 
by two factors: the background fro m solar wind (SWCX) and the 
XRB variance (Bau tz et alJl2009h . The latter was partially miti- 
gated exploiting XMM-Newton or Chandra observations. 
In the previous Section we described the procedure we devel- 
oped to quantify the impact of these two factors on Swift XRT 
temperature measurements for a sample of cluster observed in 
their central regions. Here we use the same procedures to cal- 
culate how accurate a temperature measurement performed by a 
fairly deep (300 ks) Swift XRT observation of the northern (the 
brighter) out-skirt of Abell 1795 would be. 

We used the (northern) surface brightness and temperature 
radial profile as m easured by Suzaku up to 25' and reported by 
dBautz et al.l2009 h to simulate the cluster surface brightness pro- 
file for a 300 ks Swift XRT exposure. This was done using a real 
exposure map in order to account for vignetting, CCD defects 
and different aim points and roll angles of a real observation. 
Because the Abell 1795 virial radius largely exceeds the field of 
view (r2oo=25.7'), a combination of different pointing is neces- 
sary for our purposes. We found that a mosaic of two different 
aim points with displacement of 13 ' from the cluster center op- 
timizes the observation efficiency, providing us with the max- 
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Table 4. Extraction regions and spectral analysis results of the 300 ks simulated observation. For each annulus we report the radii, 
the area, the effective exposure time (accounting for exposure maps and vignetting), the expected mean cluster surface brightness, 
the expected total events and the source/total ratio and the expected errors. We report explicitly the relative contribution of the 
statistics and the systematics terms. The latter are calculated either in the assumption that Suzaku measurement is the true value 
(HI) either considering the distribution of the possible true temperature values as expected from our simulation (H2) a s explained 
in Sect. [4] The errors are reported at 90% of confidence consistently with Suzaku published numbers (Bautz et al. 20091) 



ann. rad. 


ann. rad. 


area 


exp. 


Surf. Bright. 

erg s _I crrT 2 deg~ 2 


Tot. evt 


S/(S+B) 


T stat. error 


T syst. error (HI) 


T syst. error (H2) 




feoo] 


[deg 2 ] 


[ks] 




[%] 


[frac] 


[frac] 


[frac] 


12.6-13.3 


0.49-0.52 


1.1E-02 


8.7E+01 


2.4E-11 


1504.4 


69.72 


0.09-0.10 


0.15-0.22 


0.15-0.15 


13.3-14.0 


0.52-0.54 


1.3E-02 


8.7E+01 


2.0E-11 


1526.7 


66.05 


0.09-0.11 


0.15-0.22 


0.13-0.15 


14.0-14.9 


0.54-0.58 


1.6E-02 


8.6E+01 


1.6E-11 


1631.3 


61.37 


0.08-0.10 


0.19-0.32 


0.13-0.55 


14.9-16.0 


0.58-0.62 


1.9E-02 


9.0E+01 


1.3E-11 


1800.0 


55.64 


0.09-0.10 


0.22-0.41 


0.17-0.56 


16.0-17.4 


0.62-0.68 


1.9E-02 


1.1E+02 


1.0E-11 


1909.6 


49.41 


0.12-0.13 


0.20-0.34 


0.14-0.59 


17.4-19.0 


0.68-0.74 


1.9E-02 


1.2E+02 


1.1E-11 


2197.0 


52.07 


0.07-0.06 


0.20-0.36 


0.12-0.60 


19.0-20.6 


0.74-0.80 


1.5E-02 


1.5E+02 


9.2E-12 


2070.5 


47.39 


0.07-0.11 


0.21-0.37 


0.09-0.61 


20.6-22.4 


0.80-0.87 


1.3E-02 


1.8E+02 


9.0E-12 


2080.4 


46.78 


0.07-0.12 


0.21-0.37 


0.10-0.60 


22.4-25.0 


0.87-0.97 


1.7E-02 


1.8E+02 


6.4E-12 


2339.5 


38.39 


0.09-0.10 


0.23-0.39 


0.20-1.03 



imum exposure in the 20'-25' annulus. Starting from 1-200, we 
calculated the annuli which would contain 1000 events in the 
0.5-2.0 energy band in order to keep the statistical error at the 
level of <10% 

We used, as XRB, the same of iBautzetall d2009l) which was 
measured in the most external observed region and is very close 
and consistent with our median XRB; as NXB the NESR signal 
registered during the observation of GRB 090618 chosen as the 
most recent observation lasting more than 300 ks, 

Thus, for each annulus we calculate the total (statistical + sys- 
tematics) error in temperature measurement using the recipes de- 
scribed in Sect. [4] 

Results of this procedure are reported in Tab. |4] and shown in 
the upper panel of Fig. Q~2]in comparison with the Suzaku mea- 
surement. Systematic errors were calculated, first, assuming the 
Suzaku measurement as the true val ue of the tempera ture. This 
was done following the approach of iBautz et al.l (120091) to com- 
pare the Swift XRT expected accuracy with Suzaku (HI column 
in Tab. [4] and black error bars in Fig. [12]). The simulated spec- 
trum from the external bin is shown in Fig. [13] We found that a 
300 ks Swift-XRT observation would significantly improve the 
accuracy of the Suzaku temperature measurement in the (north- 
ern) outskirts of the cluster Abell 1795, both in terms of spatial 
binning and relative accuracy. Indeed in the annulus within 17- 
25 ' while Suzaku could measure only one single temperature 
with an accuracy of ~ 60%, we expect that the XRT observa- 
tion would be able to measure four different temperatures with 
an accuracy of ~ 40%. 

Second, we refined the error calculation, also considering the 
distribution of the possible true temperatures which could yield 
the measured value (H2 procedure, see Sect. 0J. In this case the 
upper error bars are significantly larger. This is the result of the 
not negligible bias in the high temperature (T > 5 keV) measure- 
ment in the low SB regime (Sect.|4]and FigfTTli. In this case in the 
three bins within 17-23 ' we expect an accuracy of ~ 60% simi- 
larly to the Suzaku one, while in the last bin (23-25 ') the upper 
error bar is ~ 100% of the measured value (H2 column in Tab. [4] 
and red error bars in Fig. [12]). We note that the impact of the 
bias in the high temperature (T > 5 keV) measurement in a low 
SB regime on temperature accuracy, has never been quantified 
for Suzaku. Even in the unrealistic assumption that this is com- 
pletely negligible, the Swift XRT observation would improve the 



accuracy of the temperature profile significantly narrowing the 
spatial binning up to ~ R2oo- 

In a spatially resolved spectral study of a low surface brightness 
source, a telescope reaches the limit of its capabilities, when a 
deep observation allows to map the source at the best of its an- 
gular resolution (bins >3 times the HEW, to avoid PSF mixing) 
and, at the same time, in each single bin, the collected photons 
are enough to make the statistic error negligible with respect to 
the systematics (mainly the background ones). The Suzaku tele- 
scope, in its observation of Abell 1795, almost reached its limits: 
in fact, most of the error in the temperature measurement, is due 
to background variance, while the angular resolution can be only 
slightly improved as the telescope HEW is 2'. On the other hand, 
we showed that, with a 300 ks observation, Swift-XRT could 
reach the same level of accuracy with an angular resolution of ~ 
2'. This leaves a wide range of improvement: indeed, doubling 
the observation exposure would allow to halve the spatial bin- 
ning, as the HEW is 18". 



6. Discussion and conclusion 

In the outer regions of nearby clusters, the ICM emission is only 
a small fraction of the whole signal collected by the detector. 
The regime, where the background systematics, affect the spec- 
troscopic measurements much more than the statistical error, is 
easily reached. While XMM-Newton and Chandra are not suit- 
able for this kind of observation due to the high level of par- 
ticle background, in all the published works presenting Suzaku 
observations of cluster outskirts the evaluation of the XRB and 
its varia nce is the main iss ue. Different approaches h ave been 
pursued. IBautz et al.l d2009l) and iHoshino etall d2010h . in their 
studies of Abell 1795 and Abell 1413 respectively, used the 
signal from external regions at ~ 1.2 xr2oo as XRB, in the as- 
sumption that the cluster emission is negligible at that distance 
from the center. To stud y the temperature profile of PKS0745- 
19, iGeorge et al.l ([2009) used the Lockman Hole observation 
which was performed just few days before the cluster observa- 
tion. Interestingly they found sign ificant emission from the clus- 
ter ICM at distance > 1.5 xr?nn. iKawa harada et al. (2010) used 
the two closest observations (at ~ 8° of distance) among the ones 
suitab le as blank fields in th e Suzaku archive. iReiprich et al.l 
(2009) and iFuiita et all (12008) used the classical models from 
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Fig. 12. Simulation of 300 ks observation of the northern skirt 
of Abell 1795. Upper panel: Surface brig htness of the cluste r. 
Green points are from Suzaku observation dBautz et alJ ((2009)). 
Grey points are from Swift XRT observation (already shown in 
Fig. Black points would be the result of 300 ks observation. 
Lower panel: The surface brightness profile of the cluster as 
observed by Suzaku, with the same color code. XRT spatial bins 
in the simulation are chosen in order to have a minimum of 1000 
source counts in the 0.7-2.0 keV band. 
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Fig. 13. Simulation of the 300 ks of the last annulus. 



literature to fit the cluster spectrum together with the XRB (us- 
ing a different choice of free parameters). 
In this paper we presented for the first time the analysis of the 
Swift-XRT observations of a sample of 6 galaxy clusters. We 
measured the temperature profiles as far as ~ 0.5 R200 for all (but 
Coma) sample clusters. To estimate the cosmic background we 
used a statistical approach, modeling a representative sample of 
blank fields, with a sky-coverage of ~ 15 deg 2 . We used the blank 
field median spectrum as XRB model to fit our cluster spectra. 
We calculated the systematics of this approach by simulating re- 



alistic clusters with different temperatures and surface bright- 
nesses summed to real XRB data extracted from blank fields 
with different sizes and exposure times. With this approach in the 
systematics calculation, exploiting a statistically fair sample of 
BF, we directly accounted for the XRB variance (for both CXRB 
and GXRB ). Moreover we presented a new way to calculate the 
uncertainties in the temperature measurements significantly re- 
fining the current approach in literature. Given a measure T m , we 
accounted for the possibility that T,„ is produced by a T, + T m . 
This allowed us to realistically simulate the temperature mea- 
surement in the outer regions of Abell 1795 which would be 
provided by a deep XRT observation. We showed that, thanks 
to an unprecedented combination of low background, good PSF 
the Swift XRT would be able to significantly improve the current 
accuracy of the temperature measurements in the outer regions 
of nearby clusters. 

The ideal telescope for cluster outskirts observation would be 
a large grasp (wide field and large collecting a rea) and low 
backg round telescope such the proposed WFXT dMurrav et al.l 
2010). In the next decade eRosita will be the only m ission op- 
erative with these characteristics (Pre dehl et alj 120101) . a grasp 
10 times larger than XMM (100 times larger than Swift-XRT). 
Interestingly, eRosita, with an effective area of - 1500 cm 2 at 1.5 
keV (~10 times larger than XRT) and an expected background 
of ~ 9 counts s _1 deg 2 (~10 times larger than XRT) will have 
the same source / background ratio of the Swift-XRT when ob- 
serving extended sources. If these numbers will be confirmed 
in flight, and the NXB will be reproduced with the same accu- 
racy (< 3%), at a given value of surface brightness, eRosita will 
have systematic errors on temperature measurements which will 
be very close to the ones we found for Swift-XRT. In this case, 
our proposed XRT observations would represent a pilot for the 
eRosita mission; on the other hand, if the NXB of eRosita will 
be higher than expected (eRosita will be the first X-ray telescope 
positioned in L2) or it will be impossible to reproduce it with the 
same accuracy, the XRT observation would remain the only way 
to improve our knowledge of the cluster outskirts physics at least 
for the next decade. 
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